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Using ab initio density functional total energy and molecular dynamics simulations, we study the 
effects of various forms of nitrogen post deposition anneal (PDA) on the electric properties of hafnia 
in the context of its application as a gate dielectric in field effect transistors (FET). We consider 
the atomic structure and energetics of nitrogen containing defects which can be formed during the 
PDA in various N-based ambients: N2, Nj, N, NH3, NO, N2O. We analyse the role of such defects 
in fixed charge accumulation, electron trapping and in the growth of the interface Si02 layer. We 
find that nitrogen anneal of the oxides leads to an effective immobilization of native defects such as 
oxygen vacancies and interstitial oxygen ions, which may inhibit growth of silica layer. Nitrogen in 
any form effectively incorporates into the pre-existing oxygen vacancies and, therefore may decrease 
the concentration of shallow electron traps. However, nitrogen in any form is unlikely to significantly 
reduce the fixed charge in the dielectric. 

PACS numbers: 71.55.Ht, 81.05.Je, 71. 15. Mb 



I. INTRODUCTION 

Transistor scaling, which enables continuous increase 
of performance of integrated circuits, has been generaUy 
facilitated by a reduction in the thickness of the gate 
dielectric in typical metal oxide semiconductor field ef- 
fect transistors (MOSFET). However, a further decrease 
of the thickness of conventional silicon dioxide (Si02)- 
type dielectrics leads to an unacceptable rise in the gate 
leakage current, and hence, power consumption. Several 
possible solutions are under consideration^.:, but one of 
the most attractive, which retains standard MOSFET 
design, is to replace Si02 with a material of higher di- 
electric constant (high-fc). The resulting increase in ef- 
fective capacitance means that a thicker gate dielectric 
layer can be used, reducing gate leakage current, while 
providing comparable performance to a much thinner 
Si02 layer. The choice of specific high-fc material is 
subject to rather strict electrical and process integration 
requirements^ eliminating most of the potential candi- 
dates, and presently hafnium oxide-based dielectrics are 
considered to be the most promising from a practical 
standpoint. 

In spite of recently reported successes^, the perfor- 
mance of high-fc transistors requires further improvement 
to meet the industry needs for most of the potential appli- 
cations. Most critically, high-fc devices suffer from poor 
channel mobilit^ f^r'^-^'Z and instabilities of the threshold 
potential (Vt)^- It is understood that these drawbacks 
result from the high density of structural defects in the 
high-fc films^, in particular, defects acting as electron 
trapsiiJ^. These defects were shown to form localized 
states in the band gap, whose occupancy may vary with 
the external chemical potentia]!i^4i^*i&. It was generally 
expected that a post deposition anneal (PDA) would re- 
duce the density of intrinsic defects. However, for high-fc 



films deposited on silicon, PDA usually results in a de- 
crease of the effective capacitance of the gate stack. This 
decrease was attributed to a spontaneous growth of the 
SiOa; buffer layer between high-fc and silicon substrate 
during the high-fc film deposition and subsequent anneals 
(600 - 1000 C), facilitated by the oxygen excess. The low 
quality of the resulting interfacial SiO^; layer negatively 
impacts transistor characteristics^^. In order to suppress 
the suboxide growth, nitrogen rich ambients are prefer- 
able in the PDi^. In general, nitrogen incorporation into 
the Hf-based dielectrics was shown to provide several ad- 
vantages - reducing boron penetration from the p-type 
poly silicon gate electrodes, improving electrical perfor- 
mance and increasing stability of high-fc devices. Be- 
sides N2 PDA, several other options were used to increase 
nitrogen incorporation in Hf-based oxides, including re- 
mote nitrogen plasma (N^)^^, ammonia (NH3), nitrogen 
monoxide (NO), and N2O, as well as nitrato-CVD depo- 
sition processesiSiSS. However, the lack of understand- 
ing of forms of nitrogen incorporation and its effect on 
the electronic properties of the dielectric hampers further 
progress. 

In this paper we theoretically investigate the incor- 
poration of N2, N, NO, NH3 species in different charge 
states into the ideal and defective Hf02 crystal. Specifics 
of the structural and chemical composition of real films 
depend strongly on the deposition and PDA conditions. 
In particular, a thin film deposited at relatively low tem- 
peratures (typically below 500 K) is in a metastable 
porous low density state, often characterized by substan- 
tial non-stoichiometry. Significant efforts have been re- 
cently made to resolve the atomic structure of the de- 
posited filmsSLSa. After PDA in oxygen (T > 770 K), 
hafnia films on Si adopt the monoclinic structureS^'SiiS^. 
This is in contrast with zirconia which may crystallize 
into either monoclinic or tetragonal structure depending 
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on the film thickness and stoichiometrjiSi. In this paper 
we leave the problem of polymorphism aside, and focus 
on monoclinic hafnia. 

In broad terms, an optimal post deposition anneal 
should stipulate: 

• formation of dense (preferably amorphous) oxide 
films with controlled thickness, high dielectric con- 
stant and high thermal stability; 

• reduction of fixed charge and electron traps in the 
bulk of the film and at the interface with Si (high 
electrical stability) 

• an effective control over the growth of the Si02 
buffer layer; 

• inhibition of the diffusion of charged species to- 
wards the interface and into the channel region. 

In order to assess the effectiveness of different anneals, 
we assume the following physical picture. After deposi- 
tion and PDA, a hafnia film is crystallized into the mon- 
oclinic phase which may still contain a significant con- 
centration of oxygen vacancies and interstitials. In the 
significantly non-stoichiometric films either oxygen inter- 
stitials (oxygen excess) or oxygen vacancies (oxygen de- 
pletion) will dominate. Nitrogenous species introduced 
by the PDA, diffuse into the bulk oxide and take part 
in the following (not necessarily sequential) processes: i) 
incorporation in the lattice interstitials or replacement 
of lattice oxygen ions; ii) dissociation into other species 
both in interstitial and in regular lattice sites; Hi) reac- 
tion with the existing oxygen vacancies and interstitial 
oxygen species. 

The nitrogen defects thus formed may be neutral with 
respect to the lattice or themselves form fixed charge 
or charge trap centres. This would depend on the po- 
sition of the corresponding defect levels with respect to 
the band edges of Hafnia and Silicon and with respect to 
the Fermi-level (see Fig. The latter in the n-channel 
devices is located near the silicon conduction band min- 
imum (CBM). 

All the electrically active defects in the oxide layer can 
be qualitatively classified by the position of their levels 
with respect to the silicon band edges: 

A ) Defects whose levels are resonant with the Si con- 
duction band. Under zero electric field and in the ther- 
modynamic equilibrium these levels are empty. However, 
such states are available for resonant tunneling at non- 
zero gate voltage, and thus may serve as electron traps. 

B) Defects whose levels fall into the silicon band gap. 
The electron occupancy of such defects depends on the 
position of the Fermi level. These deep defects are re- 
sponsible for the Fermi-level pinning, and thus, for the 
high threshold voltages. Also they may contribute into 
threshold potential instability and, to some extent, to a 
fixed charge problem, since the charge relaxation time on 
such defects may be macroscopically slow. 




FIG. 1: Schematic figure showing classes of defect levels with 
respect to hafnia and silicon bands. The band offset for silicon 
is taken from refi— *— . 

C) Defect states resonant with the silicon valence band 
are expected to be occupied, and, depending on their 
charge with respect to the lattice, they may become a 
major source of the fixed charge. 

Therefore, the efficiency of a specific PDA could be 
judged by the integral effect of the resultant nitrogen 
contained species on the bulk charge trap density, and on 
the mobility of oxygen, which may occur via the vacancy 
or interstitial mechanisms^^. 

Although the model outlined above is certainly simpli- 
fied, it will form a framework for systematizing the effect 
of various nitrogenous species on the properties of mon- 
oclinic hafnia. It can be then refined to include other 
important effects e.g. interface diffusion and segregation, 
impurity clustering etc. 

In the next section we give the details of the calcula- 
tions. The results of the modeling of the effect of different 
anneal gases are presented in section HI and discussion 
is given in the last section. 



II. CALCULATION PROCEDURE 

All the calculations were performed using the VASP 
code^ii^, implementing spin-polarized Density Func- 
tional Theory (DFT) and the Generalized Gradient Ap- 
proximation of Perdew and Wan^. known as GGA-II. 
The plane wave basis set was used in conjunction with 
the ultra-soft pseudopotentials of Vanderbilt type^O. The 
standard pseudopotentials for hydrogen, hafnium, ni- 
trogen, and oxygen atomic cores were generated with 
charges -|-l,-|-4,-|-5, and -1-6 respectivelji^. 

Bulk hafnia at atmospheric pressure and low tempera- 
ture has monoclinic symmetrySSiS (space group P1\jc). 
At T = 2000 K it undergoes martensitic phase trans- 
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formation into the tetragonal phase (space symmetry 
P4:2/nmc), and above T = 2870 K the cubic fluorite 
structure is most stable (space symmetry FmSm). The 
monoclinic structure is characterized by the two non- 
equivalent anion sublattices (Fig. [21). In one, oxygen 
ions arc 4-fold coordinated, and in another they are 3- 
fold coordinated. The equilibrium 0-Hf distances range 
between 2.14 and 2.24 A in the former and between 2.05 
and 2.14 in the latter sublattice. These structural pa- 
rameters are well reproduced using DFT and pseudopo- 
tentials described above and are discussed in detail in 
refii&. 

The calculations were made using a 96 atom unit cell, 
which is a 2 X 2 X 2 extension of the 12 atom mono- 
clinic primitive unit cell. We used two fc-points in the 
irreducible part of the Brillouin zone and a cutoff energy 
of 400 eV for atomic relaxation and total energy calcu- 
lations. The total energies of the charged systems were 
corrected for the spurious electrostatic interaction arising 
from the periodic boundary conditions^^i^. 

To test the stability of atomic configurations optimized 
in static calculations and to extract dynamical informa- 
tion, such as diffusion paths and vibrational frequencies, 
we carried out an extensive set of Born-Oppenheimer 
molecular dynamics (MD) simulations using VASP. In 
these simulations we used k=0 and 300 eV energy cutoff 
and a constant volume, energy and number of particles 
(NVE) algorithm with a time step of 1 fs. This assured 
a total energy drift of less than 0.05 eV at T = 600 K 
for a typical MD run of 5 ps. The supercell volume in all 
calculations was fixed to its K theoretical valuei^. 

For each charge state g of a candidate defect (D) 
we first determine its geometric structure and forma- 
tion or incorporation energy {G{D'^)) using the static 
DFT calculations. Given the macroscopically long de- 
fect diffusion and charge relaxation timesSi, we believe 
that the defect equilibrium Free energies do not necessar- 
ily reflect relative concentrations of various defects (even 
if the Fermi level and atomic chemical potentials were 
well defined). Thus, we do not discuss free energies, 
but rather assume initially that all energetically stable 
charge states are possible, and then discuss most proba- 
ble charge states considering there relative electron affini- 
ties of the defects. 

The position of defect electrical levels with respect to 
the CBM of hafnia and silicon determines the stability 
of a defect in a particular charge state as well as its 
ability to trap or release an electron. These properties 
are quantified by the vertical and relaxed defect ioniza- 
tion energies (/p(I?^)), and vertical electron {xe{D'')) and 
hole (x/i (£*')) affinities, which are calculated using the 
approach described in ref.^^ (see Appendix for details). 
This approach is based on standard definitions of elec- 
tron affinities and ionization energies, it employs combi- 
nations of total energies of the system in different charge 
states and conserves the number of electrons. In order 
to discuss the electron and hole trapping by defects from 
silicon, we neglect the band bending at the interface and 



the effect of applied voltage and use the experimental 
band offsets between bulk hafnia and siliconS2i2&. 

Towards the full picture of possible products of incor- 
poration of nitrogen, we examine all possible dissociation 
channels for nitrogenous molecular species in hafnia ma- 
trix. This is done by calculating their dissociation energy 
into various products. The latter is found as the differ- 
ence of total energies for the system with the associated 
defect and the joint energy of the infinitely separated 
dissociation products, corrected for double counting of 
the perfect supercell term. In further discussion we treat 
a molecular defect as thermodynamically unstable if its 
dissociation is exothermic. 



III. RESULTS OF CALCULATIONS 

As discussed previously, nitrogen-based PDA can be 
performed using a variety of gas sources (ambients), and 
depending on the chemical composition, pressure and 
temperature, different products may occur in the bulk 
film. Here we consider each PDA ambient in turn, and 
predict the likely resultant defects and their effects on 
the properties of the oxide. 



A. Molecular Nitrogen 

1. N2 incorporation 

The experimental dissociation energy of an N2 
molecule in the gas phase is ~9.9 eVM, (10.0 eV in our 
DFT calculations). Thus, no N2 dissociation is likely to 
occur at the surface, and molecular nitrogen is expected 
to diffuse inside the bulk oxide. Indeed, as we show be- 
low, even trapping of extra electrons or interaction with 
oxygen vacancies does not lead to a spontaneous dissoci- 
ation of N2 molecular species inside the oxide. 

The equilibrium configuration of the neutral interstitial 
nitrogen molecule is shown in Fig. 12 The incorporation 
energy for this configuration with respect to N2 in the 
gas phase is ~3.0 eV per nitrogen atom. It is seen in 
figures 121 and Ufa) that the interstitial N2 does not form 
a bond with the lattice oxygen, but rather incorporates in 
the interstitial space near to the three-fold coordinated 
oxygen site. This is in contrast with the O2 molecule, 
which makes a covalent bond with the lattice anionic. 

The interstitial N2 molecule has a positive affinity to 
two electrons (defect electron affinities are given in Table 
U). The localization of the first extra electron on the N2 
interstitial (N2^ radical) results in the outward relaxation 
of the nearest neighbor oxygen ions (cf. 0-N distance of 
1.96 A in and ~2.20 A in N^). The corresponding 
relaxation energy is ~1.2 eV. The second extra electron 
is also strongly localized on the molecular ion. This 
further increases the 0-N distance to 2.34 A with a relax- 
ation energy of 0.6 eV with respect to the geometry. 
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FIG. 2: Structure of neutral nitrogen interstitial molecule in 
hafnia. The ions are colour coded as follows: white - Hf, grey 
- O, black - N. 

Interestingly, a hole (an absence of the electron) in Hf02 
will not localize on the N2 interstitial, owing to the fact 
that the last occupied orbital of the molecule is well be- 
low the valence band maximum (VBM). 

TABLE L Vertical ionizational potential I{D), relaxed elec- 
tron Xe{D) and vertical hole Xh{D) affinities (in eV) of defects 
in different charge states in hafnia. The subscript O indicates 
the substitution in the 3-coordinated oxygen site. 



D 


IpiD) 


Xe(D) 




N" 


5.0 




0.7 


N° 


5.3 


4.2 


0.3 


N+ 


5.6 


4.1 


0.1 






5.0 




Nr 


4.2 




1.5 




4.3 


3.4 


1.3 


N§ 




3.3 




NO" 


4.8 


4.4 


0.7 




5.2 




0.5 




5.6 


4.6 


0.1 


K 




5.2 




(N2)o 


3.9 




1.7 


(N2)S 


4.4 


3.3 


1.1 


(N2)S 




3.8 





2. Thermodynamic considerations on the likely charge 
states 

We have established that interstitial nitrogen 
molecules may have multiple charge states, whose 
occurance at thermodynamic equilibrium is defined 
explicitely by the electron chemical potential. The latter 
depends not only on the external potential, but also 
on the concentrations and electrical levels of all other 
defects present in the lattice. In this case, the problem 
of the possible defect charge states may be evaluated by 
considering relative electron affinities of various defects, 
and answering the question: given the fixed number of 
electrons, on which of the defects they are more likely 
to reside? This question can be resolved by considering 




FIG. 3: Charge density 2D slices around various nitrogen 
defects in hafnia: (a) plane through an interstitial nitrogen 
molecule and a three coordinated oxygen atom; (b) plane 
through a neutral nitrogen interstitial, three coordinated oxy- 
gen and hafnium; (c) plane through a single positively charged 
nitrogen interstitial and two lattice oxygens. 
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TABLE II: Electron exchange reactions for defects in m-Hf02. 
The energies are calculated in the assumption of infinitely 
separated defects. _E > indicates exothermic reaction along 



the arrow. 


i\o 


Reaction 


Energy (eV) 


1 
1 




1 1 
i.i 


o 
z 




A 

Z.4 


Q 

o 


7\rO 1 T/0 AT-— 1 T/" + 

iV2 + V ^ + V 


u.o 


4 


7\rO 1 T/{) _v Tvr 1 \r-\ — V 


A 7 
U. / 





iV2 + V ^ iV2 + \/ 


A 1 
U.i 


D 


iV2 + V => I\2 + V 


A a 
U.D 


I 






o 

O 


AT-Q 1 T/0 1. AT— 1 T/ + 


1 /I 

1.4 


Q 


A/" n"l~ _i_ T/"^ ^ AT /T^ _L + 
i\ ~r V =p i\ 1I2 + V 


i.O 




0"+ 1 T/0 _v rr— 1 T/H — h 


A A 

-U.4 


11 


AT— 1 /"\ — , AT 1 /~iO 

A/2 + u => A/2 + u 


-0.5 


12 


TV/-— 1 — , TvrO 1 /~i 

JM2 ~\~ U => A/2 + c 


1.2 


13 


A! 1 > A rO 1 

A'2 + A'2 + 


1.7 


14 


A* + A + C 


-0.2 


15 


a" + A" + 0° 


0.3 


16 


N- +0- ^ N° + 0~- 


0.3 


17 


N++ +0— ^ N+ +0- 


0.6 


18 


N++ +0- ^ N+ +0° 


1.2 


19 


2V+ ^V" + V++ 


0.3 


20 


20- => 0° + 0— 


0.6 


21 


2N^ ^ + N^' 


0.1 


22 


2N° ^ A+ + A~ 


0.1 


23 


2NHi ^ NH+ + NH^ 


2.8 


24 


2H° ^ H+ + H- 


2.7 



the electron exchange reactions between the pairs of 
infinitely separated defects as presented in Table HTl 

Consider for example the oxygen Frenkel pair (an in- 
terstitial atom and the vacancy). It is evident from reac- 
tions 1^1,2) that in the overall neutral m-Hf02 the zero 
temperature equihbrium state corresponds to the dou- 
bly positively charged vacancies and doubly negatively 
charged interstitial oxygen atoms. Next, consider N2 in- 
terstitial molecules. Reactions ((^3-6) suggest that they 
are more electronegative then anion vacancies. There- 
fore, in oxygen deficient films (high concentration of va- 
cancies) nitrogen molecules will be predominantly doubly 
negative, while the vacancies remain doubly positive. 

In the case of oxygen excess (large concentrations of 
interstitials) , the relative charge state is defined by the 
energy balance in the reactions 1^11-13). It follows that 
interstitial oxygen ions are more electronegative than ni- 
trogen molecules. However, as discussed above, is 
not stable, so nitrogen interstitial molecules will remain 
largely neutral in this case. 

Finally, reaction (jn|21) suggests that neutral and dou- 
bly negative N2 interstitials are only marginally prefer- 
able over the N^. 

The energetics of the electron exchange reactions con- 
sidered in Tab. ^ allows to narrow down the range of 
relevant association/dissociation reactions between vari- 
ous pairs of defects considered below. 



TABLE III: Energetics of the association reactions of nitrogen 
and the interstitial oxygen in m-Hf02. E > denotes an 
exothermic reaction along the arrow. 



INO 


Reaction 


— — 

Energy (eV) 


1 


ZU U2 


1 A 

-l.U 


2 


o° + o- ^ O2 


-0.3 


3 


20- ^O^' 


0.7 


4 


o° + o-- ^02~ 


0.0 


5 


(A2) + 0" => (A2O) 


-0.9 


6 


{N2)+0° => {NO)+N 


-2.1 


7 


N" + 0" ^ (NO) 


1.0 


8 


A" +0° ^ {Noy 


1.2 


9 


N'' + 0- ^ {Noy 


1.6 


10 


A+ + 0° ^ {NO)+ 


1.5 



3. N2 reactions with interstitial oxygen species 



First, we observe that the interstitial oxygen molecules 
(O2) are either marginally stable or unstable in any 
charge state (Table 1-4). Therefore, we conclude that 
excess oxygen is present only in the form of atomic inter- 
stitials. As demonstrated in refsii^i2L22i2S, an interstitial 
oxygen atom forms a dumbell type configuration in silica, 
zircon, zirconia and hafnia where an oxygen atom effec- 
tively incorporates into the Me-O-Me bond (Me is Si, Zr 
or Hf ion) forming a Me-O-O-Me configuration with a 
0-0 covalent bond (1.5 A in Hf02). 

Second, it follows from the calculations that both N2 
and oxygen interstitials may accept charge states 0,-1,-2 
with respect to the lattice (Table IJ). Their likely as- 
sociation product, the N2O interstitial, is stable only 
in charge states (+2, +1,0). Next, it follows from the 
electron exchange reactions shown in TablelTll 11-13 that 
oxygen interstitials are more electronegative than nitro- 
gen molecules. Therefore, relevant reactions involving N2 
and O must have nitrogen molecules in the charge state 
not more negative than that of the oxygen. Disregard- 
ing also reactions involving unstable initial or final charge 
states leaves only one possible association reaction (Table 
inn s'), which is endothermic. Therefore, the formation of 
neutral N2O interstitials is unfavorable. Similarly, it is 
kinetically unlikely that interstitial NO is formed from 
N2 molecules, since all reactions involving N2 dissocia- 
tion (e.g. IIIII6) are strongly endothermic. The fact that 
the N2 molecules do not bond to the interstitial oxygen 
atoms suggests that molecular nitrogen is ineffective in 
inhibiting diffusion of highly mobile oxygen interstitials^^ 
and in preventing them from migration towards the in- 
terface. This finding is consistent with experimental ev- 
idence that oxygen abundant deposition techniques lead 
to the formation of the Si02 buffer layer at the IIf02/Si 
interface, and that the N2 PDA does not significantly 
affect its thicknessiiSi. 
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TABLE IV: Energetics of the reactions of some interstitials 
with a neutral 3-fold coordinated oxygen vacancy. The sub- 
script O denotes the specie incorporated into the oxygen va- 
cancy. Energies are calculated in the assumption of the in- 
finitely separated vacancy and interstitial in the initial (left 
hand side) state. E > Q denotes an exothermic reaction along 
the arrow. 



No 


Reaction 


Energy 


1 




5.1 


2 


N- + V++ ^ {N2)+ 


4.7 


3 


N^- + v+^ mo 


5.4 


4 


N- + V++ => N+ 


2.2 


5 


N- +V+ ^ iVg, 


4.3 


6 


N° +V+ ^ N+ 


3.3 


7 


N° + V° ^ 


5.7 


8 


N- ^ N- 


6.1 


9 


{NH2)+ + l^" ^ {NH2)+ 


6.9 


10 




6.6 


11 


iH)+ + V''^H+ 


2.7 



4.- N2 reactions with oxygen vacancies 

In the case of substantial oxygen vacancies concen- 
tration in the film, the N2 interstitial species can in- 
teract with the vacant oxygen sites and passivate them. 
Based on the previous calculationsi^ii&, we assume three- 
coordinated oxygen vacancies to dominate. We con- 
sider the stability of nitrogenous species in the three- 
coordinated anion lattice sites with respect to the in- 
finitely separated interstitial species and the oxygen va- 
cancy. Furthermore, we recall that nitrogen interstitial 
molecule is more electronegative than the oxygen vacancy 
(see Table and discuss reactions involving doubly 
charged vacancies {V'^^) as the most thermodynamically 
probable. The corresponding reaction energies are sum- 
marized in Table IWI where the subscript O denotes the 
specie incorporated into the oxygen vacancy. Note, that 
the superscript denoting the charge state corresponds to 
the charge with respect to the lattice^ which differs from 
the molecular charge in case of the substitute molecules. 
For example, the charge localized at the N2 molecule (nu- 
clear plus electron) in the rhs of the reaction (|IVI 1) is -2, 
so this defect is isovalent to the host O ion, and thus, 
is neutral with respect to the lattice. Similarly (^^2)0 
the rhs in (|1VI 2) denotes the substitution, which is 
in a charge state -1-1 with respect to the lattice. 

One can see that most of the stable interstitial species 
react with the pre-existed anion vacancies with a very 
substantial energy gain. Reactions IIVI 1.2 correspond to 
a formation of N2 molecule at the oxygen site, which is 
structurally similar to the 0-N complex formed by the 
nitrogen atomic interstitial (see below and Fig. |SIb)). It 
is interesting to note, that the (N2)o is stable in three 
charge states -1-1,0,-1. It is evident that formation of the 
negative (N2)o substitution involves either pre-existing 
neutral vacancies (vacancies with 2 trapped electrons), 
or an electron(s) trapping by the neutral substitution. 



Although such processes are even more energetically fa- 
vorable (e.g. I1VI 3). they might be kinetically inefficient 
due to the much lower mobility of the V*' and V+ as 
compared to V+"''. 

Adding/removing an electron to/from the (N2)o 
just increases/decreases the N-N bond length and de- 
creases/increases the nearest neighbor Hf-N bonds. The 
electron affinity of N2 at the oxygen vacancy is similar to 
that of the interstitial molecule (Table P). 

The energy gain in the interstitial-vacancy reactions 
shown in Table IIVI may be compared with the forma- 
tion energy of the oxygen Frenkel pair in different charge 
states (8.1, 7.0 and 5.6 eV for the (V° - 0°), (V+ - O"), 
and (V^"*" - 0~— ) pairs respectively^^). It follows that, 
although the substitution of site oxygen ions by nitrogen 
molecules is endothermic, the substitution energies can 
be as low as 0.4 eV (as for Nj'^ molecular ions), so the 
processes involving oxygen replacement by nitrogen are 
possible at elevated temperatures. 

We note that the reactions involving different charge 
states of the vacancy can be readily calculated but they 
are likely to be even more exothermic due to a less stable 
left hand side (e.g. reaction (|IVI3). 



B. Atomic Nitrogen 

As already mentioned, a dissociation of the N2 
molecule in the gas phase costs about 10 eV. However, 
this energy is greatly reduced in the bulk of hafnia crys- 
tal. Table [V] summarizes all possible reactions which 
would lead to the dissociation of N2. In reactions (|Vl l.2) 
we see that dissociation of the neutral molecule within the 
crystal is still endothermic, requiring about 3 eV for both 
neutral and charged products. The dissociation energy 
of the Nj" molecular ion, although smaller, is still ~1.6 
eV ((33,4). Thus, even without considering the reaction 
barriers, the dissociation of the nitrogen molecule seems 
highly unlikely even at elevated temperatures. We also 
find that, similarly to the N2 interstitial species, the N2 
species trapped by oxygen vacancies are stable with re- 
spect to dissociation. The corresponding reactions l(Vl 5- 
11) refer to N2 species in the oxygen vacancy on the left 
and to one of the product species in the vacancy and 
another in the interstitial position on the right. As one 
can see, all dissociation reactions considered in Tab. |V] 
are endothermic, hence, trapping of N2 in the oxygen 
vacancy does not significantly reduce its dissociation en- 
ergy. From this, one may infer that after N2 PDA, the 
concentration of atomic or ionic nitrogen in the bulk film 
will be negligible. 

However, atomic nitrogen can be introduced into the 
system by the plasma assisted nitridationi^ (Nj), or us- 
ing metal- nitride precursors in the film depositioniS. 
Therefore, we consider next the energetics of the atomic 
nitrogen species. 
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TABLE V: Dissociation reactions and associated energies for 
N2 in hafnia. The energies are calculated in the assumption of 
the infinitely separated products. Negative energies indicate 
endothermic reaction in arrow direction. 



No. 


Reaction 


Energy (eV) 


1 


„ , 


-3.1 


I 


i\rO i ?\r+ 1 \J — 

1^2 i V + iV 


Q 1 

-6.1 


3 


^ AT" + iV- 


-2.2 


4 


Nj- ^ N- +N- 


-1.1 


5 
6 


{N2Z ^ {N-)o + N+ 
{N2t^{N")^+N'> 


-2.9 
-3.4 


7 


{N2fo^ {N+)^ + N- 


-4.4 


8 


{N2)o^{N-)^+N^' 


-2.1 


9 


{N2)o^{N°)^ + N- 


-2.5 


10 


(iV2)+^ (A+)^+A0 


-4.7 


11 




-3.6 



1. N atom incorporation 

The interstitial nitrogen atom in its ground state con- 
figuration makes a covalent bond with the three-fold co- 
ordinated lattice oxygen. This 'dumbell' structure is sim- 
ilar to that of the interstitial oxygen atom discussed pre- 
viously in zircon'^®, zirconia'^^ and hafnia^^. The struc- 
ture of the defect is shown in Fig. ^ and the correspond- 
ing electron density is depicted in Fig. Otb). The elec- 
tronic ground state is a doublet. This is in contrast with 
atomic nitrogen in vacuum, whose quadruplet ground 
state is almost 3 eV lower than the doublet. An addi- 
tional electron also fully localizes on the 0-N bond with 
an affinity of 4.3 eV, increasing the bond length by about 
0.1 A. The relaxation energy after the electron trapping 
is about 1.2 eV, compared to over 2 eV for interstitial 
oxygen. In contrast to the N2 molecule and an atomic 
oxygen, the atomic nitrogen interstitial has no affinity for 
a second electron. However, atomic nitrogen is also sta- 
ble as a positive ion. The local relaxation in this case is 
distinctly different (Fig. |3^c)). The nitrogen-oxygen pair 
is now very electron deficient, and in fact nitrogen forms 
a weak bond with a second lattice oxygen (N-0 bonds 
are 1.46 and 1.59 A). This is accompanied by an energy 
gain of ~1.2 eV. Removal of another electron forces the 
nitrogen to a symmetric position, with two strong equiv- 
alent bonds (1.40 A) to lattice oxygen sites and an energy 
gain of ~0.6 eV. Thus interstitial atomic nitrogen may 
exist in monoclinic hafnia in four stable charge states 
(+2,+l,0,-l). 



2. N reactions with interstitial oxygen species 

In contrast with the N2 , an atomic nitrogen binds to an 
oxygen interstitial with an energy gain of around 1.0-1.5 
eV depending on the charge state (Table ^^7-10). This 
process is associated with the breaking of the 0-0 bond 




FIG. 4: Structure of neutral nitrogen interstitial in hafnia. 

of the 'dumbeir and formation of the NO2 quasi-molecule 
(Fig. lib). The 0-0 distance is now ^2.2 A (cf. 1.5 A 
in the neutral dumbell^-), while the two N-0 bonds are 
^1.3-1.4 A. Interestingly, in this configuration there is 
no clear distinction between the interstitial and the site 
oxygen - the local symmetry for both oxygen ions remains 
similar. The energy gain in reactions of charged oxygen 
and nitrogen interstitial atoms is even greater (reactions 
IIIII 9.10). Note however that in contrast with N2, the 
electronegativity of the atomic interstitial is very close 
to that of the oxygen (^14-18), so N and O interstitials 
are expected to have similar charges. 

The stability of the interstitial NO molecular species 
suggests that atomic nitrogen may play an important role 
in immobilizing fast diffusing oxygen interstitials, espe- 
cially O^, thereby inhibiting the growth of the unwanted 
Si02 interfacial layer during anneal. 



3. N reaction with oxygen vacancies 

In considering reactions with vacancies we recall that 
the interstitial atomic nitrogen may accept the charge 
states +2, +1, 0, -1, the anion vacancy - +2, +1, 0, and 
the atomic nitrogen substitution: -1-1,0,-1. We also take 
into account the result that the nitrogen interstitial is 
substantially more electronegative than the vacancy, so 
the relevant reactions involve nitrogen in a charge state 
not more positive than that of the vacancy. These con- 
siderations leave five relevant reactions shown in Table 
I1VI 4-9. Similarly to N2, atomic nitrogen passivates the 
neutral vacancy with an excess energy in the range of 
2-6 eV. One should also expect the vacancy passivation 
by atomic nitrogen to be more effective kinetically due 
to a higher mobility of atomic versus molecular nitrogen. 
Notice, however, that for the same reason reactions in- 
volving V++ will be faster than those with V+ or Vq as 
discussed before. 

The nitrogen atom incorporation into the vacancy 
causes marginal outward relaxation of the three nearest 
Hf ions. A nitrogen atom at a vacancy has the highest 
electron affinity of any of nitrogen defects (see Table 0. 
The addition of an electron (in other words N~ incorpo- 
ration) causes the N-Hf bond lengths to become smaller 




than the 0-Hf distance in the perfect lattice. The cor- 
responding relaxation energy is about 0.6 eV. For N+ 
substitution the reverse relaxation occurs, i.e. the bond 
lengths increase, but with a similar energy gain. 

An interesting and potentially important feature of 
both the nitrogen atom and molecule incorporated in an 
oxygen site is that both species can exist in the negative 
charge state. This means that the oxygen vacancy with 
N or N2 in it can effectively accommodate three electrons 
and the third electron is strongly bound see (Table QJ. 

On the other hand, our DFT calculations do not pre- 
dict any electron affinity for the neutral oxygen vacancy 
in hafnia (see refs.^'*''*° for discussion). Thus the forma- 
tion of negatively charged centers is stipulated by the 
large electron affinity of nitrogen species. Together with 
interstitial nitrogen species, these centers can serve as 
deep electron traps and be responsible for the negative 
oxide charging. 

However, the existence of negatively charged nitrogen 
species may have also some positive effect as they interact 
with protons. For example, reaction I VII 5 clearly shows 
that it is energetically favorable for the interstitial N~ 
to trap the proton and form a neutral NH interstitial. 
Since atomic nitrogen has a high electron affinity, and is 
likely to exist in the negative charge state, this provides 
a method for removing excess charge in the oxide due to 
protons - a problem often seen in experimental studiesi^. 
However, the small energy gain of -0.6 eV for the reaction 
implies the possibility of some NH dissociation at higher 
temperatures. 



C. Properties of NO and N2O in hafnia 

Our calculations of the interaction of N and N2 with 
excess oxygen in the oxide have established that an NO 
interstitial binds with lattice oxygen sites. Reaction 
IIHI 5 viewed from right to left clearly suggests that N2O 
molecules are likely to dissociate into neutral or charged 
N2 species and oxygen interstitials. It is likely that nitro- 
gen dioxide (NO2) will also dissociate to produce NO and 
oxygen interstitials. Therefore, PDA using nitric oxide. 



TABLE VI: Energetics of the dissociation reactions for am- 
monia and its products in m-Hf02. The stable configuration 
of a proton corresponds to a bonding to the three-coordinated 
oxygen. E > denotes an exothermic reaction along the ar- 
row. 



No 


Reaction 


Energy (eV) 


1 


{NH3) ^ {NH2)- +H+ 


-0.8 


2 


{NHz) ^ {NH) + H2 


-1.6 


3 


{NH2)- ^ {NHy + H- 


-1.7 


4 


{NH2)+ ^ {NH)° + H+ 


-0.3 


5 


{NH) ^ N- +H+ 


-0.8 


6 


{NH)+ ^ {N)° + H+ 


0.1 



dinitric oxide, or nitrogen dioxide ambients will not be 
efficient in trapping oxygen and in preventing growth of 
SiOx buffer layer at the interface. 

According to our calculations, the NO interstitial is 
a deep centre with an ionisation potential and elec- 
tron/hole affinity similar to those of the N interstitial (c.f. 
Tabled). Being an open shell molecule, it is also stable 
in three charge states (+1,0,-1). The positive or negative 
charging of this molecule is associated with a large and 
not trivial lattice relaxation (Figure Isj. As discussed in 
the previous section, an NO interstitial bonds to the 3- 
coordinated lattice oxygen, forming an NO2 -like config- 
uration at the site. The orientation of this configuration 
is strongly coupled to its charge state. The ONO angle 
also changes from 115° for a positive molecule to 108° for 
the neutral, to 103° for the negative molecule. The relax- 
ation energy upon ionization is 1.3 eV, and that upon an 
electron trap is ~1 eV, confirming the strong structural 
rearrangement . 

NO and N2O molecules can interact with anion va- 
cancies by incorporating with their oxygen end into the 
vacancy. The structure of the resulting defects is equiva- 
lent to that of atomic and molecular nitrogen interstitials 
discussed earlier in the paper. 

Nitrogen interstitial atoms can trap interstitial oxygen 
according to reactions III1I 7-10. However, the contribu- 
tion of this sequence of reactions is likely to be small 
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since it requires an initial vacancy as a catalyst. 



D. Incorporation, diffusion and dissociation of NH3 

Apart from its usage for PDA of hafnia films, incorpo- 
ration of ammonia into oxides is an interesting example 
of a complex molecular defect. We shall start with the as- 
sumption that ammonia incorporates into the monoclinic 
hafnia lattice intact. However, we shall see, that such 
states are metastable and more stable dissociated config- 
urations exist. The fact that the molecule is hydrogen 
rich prompts the possibility that some hydrogen-related 
defects can be formed. In fact protons have been impli- 
cated in hole trapping in hafnia films and at interfacesi^. 
Therefore we pay particular attention to possible disso- 
ciation paths of ammonia in hafnia. 

In contrast to the molecular species discussed previ- 
ously, ammonia is a nucleophilic molecule and it binds to 
lattice Hf ions. The stable configuration corresponds to 
the ionic Hf-N bond aligned approximately along (011)- 
type directions (Fig. IH^a)), with the C3 axis of the am- 
monia molecule aligned with this bond. The equilibrium 
Hf-N bond length is approximately 2.1-2.2 A . The next 
nearest neighbor Hf(2) ion lies approximately on the same 
axis and is separated from the nitrogen core by only 2.5- 
2.8 A . 

To understand the mechanism of diffusion of ammonia 
inside the lattice one needs to analyze its vibrations. The 
molecule's thermal motion involves primarily four dis- 
tinct groups of modes: stretching of the Hf(l)-N bond 
(S), libration of the Hf(l)-N bond (L), rotation of the 
NH3 around its C3 axis (R) , and the ammonia's so-called 
'umbrella' mode (U) (Figure El^a)). Two nearly degen- 
erate L-type modes have relatively low frequencies and 
manifest themselves by large amplitude 'swings' of the 
NH3 molecule, in which C3 axis of the molecule remains 
coUinear with the Hf(l)-N bond. These 'swings' are cou- 
pled to the R-type rotations of similar frequencies, so 
the directions of NH bonds of the molecule may adjust 
to a local environment. This L-R type thermal motion 
remains reversible until the C3 axis is aligned with the 
Hf(l)-Hf(2) direction. Once the Hf(l)-N-Hf(2) bonds 
are aligned, the S-mode softens, so a swap between the 
short and long Hf-N bonds becomes possible. As seen 
in the inset of Figure EJa), such molecular hops are con- 
trolled by S and U vibrations of ammonia and involve 
only a small nitrogen displacement of under 0.4 A. The 
activation energy of the umbrella mode is ~0.1 eV (gas 
phase datai^), so frequent NH3 hops between few shal- 
low minima are readely observed in the MD simulations 
even at T=300K. An alternative mechanism of ammonia 
hops has transient states between the oxygen bridged Hf 
atoms (Hf(l) and Hf(3) in Fig. IHl^a)). This involves direct 
S-L type couphng (see inset Fig. I^a)) and has a higher 
activation energy. Despite being fairly frequent (at least 
one event in 1 ps at T=3GGK), S-U and S-L type hops do 
not contribute to molecular diffusion, since in such hops 



the molecule does not leave the original Hf02 interstitial 
cage. 

We have been unable to find a diffusion path for the 
whole NH3 molecule. In fact, further calculations show 
that ammonia molecule in the m-Hf02 (Fig. El^a)) is 
metastable and dissociates into NH^ and a proton. The 
compact (NH^-H+) pair is best described as an NHj" ion 
replacing an oxygen ion in a 3-coordinated site, while the 
displaced oxygen traps the proton and relaxes towards 
the volume interstitial position (Fig. Efb)). The energy 
of this configuration is by ~ 0.7 eV lower than that of 
the associated molecule. However, the energy barrier be- 
tween the two configurations is substantially higher than 
thermal energies. This is reflected by the fact that the 
molecular NH3 configuration at 300K is stable during the 
entire MD run (at least 6 ps). However, at T = 600K, 
which is still below usual PDA temperatures, the disso- 
ciation occurs typically within 1-2 ps. 

Stability of the dissociated configuration (Fig. Elb)) 
is stipulated by the fact that the OH" (or NHj") bonds 
are stretched by the large electrostatic potential gradi- 
ent near the anion site. Dipole fields of the NH2^ and 
0H~ shearing a single site, mutually counterballance the 
crystal field effect, and stipulate stronger proton bonds 
and local neutrality. Consequently, although the initial 
NH^-H+ dissociation between states a and h is exother- 
mic, an infinite separation of H+ and (NH)^ is endother- 
mic by 0.7 eV ( Table IVTI l'l. The latter can be achieved 
e.g. by proton diffusion away by hopping between the 
3-coordinated oxygen ions (Fig. Et^b)). An alternative 
dissociation path, resulting in the formation of a neutral 
molecular hydrogen and NH fTable rvTl 2). is substantially 
less energetically favorable. 

When separated from the proton, the NH^ radical may 
change its charge state when an external potential is ap- 
plied. Our calculations suggest that this molecule can 
be either positive or negative, whilst the neutral config- 
uration is thermodynamically unstable at any value of 
Ej;^ (Table ini 23). The same apphes to hydrogen ||ni24), 
which is only stable as a proton or a negative ion in hafnia 
(see referencea^Mi for discussion on hydrogen) . Depend- 
ing on the charge state, the dissociation of NH2 may lead 
to the formation of a neutral nitrohydratc and a negative 
or positive hydrogen ion fTable lVll 3.4). Both reactions 
are endothermic by 1.7 and 0.3 eV, respectively. Despite 
the apparently low dissociation energy, our MD simula- 
tions at T = 500 K shows that the NHj ion is stable - 
indicating that the dissociation barrier may be high. 

Depending on its charge state, the NH2 molecule binds 
to the lattice quite differently. The NH^ ion is iso- 
electronic to the 0H~ radical. Therefore, as expected, 
its most stable configuration is similar to the interstitial 
OH" centre, as reported for monoclinic zirconi£ii&. NH^ 
binds to three Hf ions with a bond length of around 2.3 
A (Fig. Efa)). In contrast, the NH^ ion makes an ion- 
covalent bond with a 3-coordinated oxygen ion (0-N dis- 
tance 1.4 A ) and an ionic bond with only one of the 
cations (N-Hf distance 2.2 A) (Fig. ^h)). 
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FIG. 6: Equilibrium configurations for the ammonia molecule a) intact NH3 b) locally dissociated form: NH^ + H+ (more 
stable by 0.7 eV). The oxygen ion accepting the proton (marked as I) moves off the lattice site indicated by 'X'. The insets 
of the figure a) illustrate the vibrational modes involved in 'umbrella' mode assisted jump (left) and 'swinging' mode assisted 
jump (right) of the NH3 molecule. 




FIG. 7: Equilibrium configurations of the NH2 interstitial (a) 
NH^; b) NH+. 



Further dissociation of NH2 may lead to the formation 
of NH species. Calculations predict NH" to be stable, 
while NH~ is unstable with respect to electron donation 
to the CB, and NH+ is unstable with respect to proton 
release (TableEl6). NH° molecule forms a dumbell con- 
figuration with a 3-coordinated lattice oxygen similar to 
N interstitial (N-0 distance is 1.47 A ). The equilibrium 
angle ONH is 107°, but is soft. Finally, dissociation of 
NH° is energetically unfavorable, as discussed previously 
(Table Ell 5). 

The presented results suggest that the likely final prod- 
ucts of ammonia PDA are NH^ , NH*^ as well as hydrogen. 
The latter is stable either as a proton or as a negative ion. 
The detailed dynamical behavior of hydrogen will be dis- 
cussed elsewhere. NH^ and NH" species are unable to 
trap interstitial oxygen. This agrees with experimental 
studies showing an increase in interface Si02 thickness 
during NH3 PDA. In addition, NH^ and NH" (but not 



NH^ ) can react exothermally with the existing anion va- 
cancies (see Table llVl 9. 10). The resultant defects create 
electronic levels near or inside the valence band of hafnia. 
Therefore, unlike atomic and molecular nitrogen species 
trapped in oxygen vacancy, NH^ and NH" substitutions 
are stable in only one charge state. 

IV. DISCUSSION 

Our calculations identified the most stable nitrogen 
containing species and the corresponding energy levels 
with respect to the band edges of the monoclinic hafnia. 
It is revealed that most of the stable nitrogen contained 
products (as well as intrinsic defects) may have multiple 
charge states depending on the availability of the elec- 
trons in the system. In the MOS architecture the main 
source of electrons in the dielectric is the channel silicon 
and poly-Si or metallic gate. Therefore, to identify the 
most likely charge states, the calculated electric levels 
of defects must be aligned with the silicon bands. Such 
an alignment is illustrated in Figure |H1 where the elec- 
tric levels are adjusted to the experimental band gap en- 
ergy as described in the Appendix, and the experimental 
band off-sets are used for the valence and the conduction 
bands2Si2&. Here we consider main effects caused by these 
defects. 



A. Structural effects 

Our calculations predict that all the neutral and nega- 
tively charged stable nitrogen species (N2, N, NH2" , NH", 
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and H ) incorporate into oxygen vacancies with a sub- 
stantial energy gain. However, this process is diffusion 
limited. Since the diffusivity of N2, NH2, NH°, and H" 
is unlikely to be high, the process of vacancy passivation 
by molecular nitrogen or ammonia may be incomplete, 
leaving significantly non-uniform concentrations of inter- 
stitials and anion vacancies. In contrast, atomic nitrogen 
is expected to be rather more mobile resulting in more 
effective passivation of oxygen vacancies. A similar argu- 
ment applies to the processes involving anion vacancies 
in charge states -1-1 and 0. They are significantly less 
mobile than the bare V''""'" due to the electrons localized 
inside. Consequently, the time to attain thermodynamic 
equilibrium might be long, and the subtitution by the 
PDA agent in practical anneals may be incomplete. 



B. Electron trap effects 

Our calculations put the oxygen vacancy electrical lev- 
els into the silicon band gap (Figure |H1) , which would 
make them type B defects. However, this attribution 
is a result of our calculation procedure in which the 
DFT band gap correction is applied to electron affini- 
ties of all defects irrespectively of their origin (see Ap- 
pendix). Such approximation may become inadequate 
for the doubly charged oxygen vacancy whose (unoccu- 
pied) electronic state splits down from the CBM due to a 
strong outward displacement of the three nearest neigh- 
bour hafnium ions. Thus, a vacancy state originates from 
a hafnia conduction band, and its electron affinity ought 
to be calculated accordingly. Assuming that the vacancy 
level must be shifted up rigidly with CBM, one would 
obtain relaxed and unrelaxed electrical levels V++/V+ 
at 0.5 eV and 1.3 eV below CBM. These values give a 
lower limit for the V"'""*' vertical electron affinity. For the 
sake of consistency we retain the values obtained by the 
defect independent procedure outlined in the Appendix, 
but acknowledge that anion vacancies may belong to a 
type A, i.e be resonant with the Si CB. Therefore, they 
may act as shallow traps. Accepting this as a working 
assumption, the role of nitrogen PDA with regard to the 
shallow trap problem reduces to the effectiveness of the 
incorporation of nitrogen species into anion vacancies. 

Nitrogen incorporation into the oxygen vacancy sites 
leads to the lowering of the defect levels, making them 
less active as electron traps. In this respect, all the con- 
sidered PDA will reduce the shallow trap concentration, 
although to a different extent as discussed before. 

However, another shallow trap candidate could be a 
small radius electron polaron in the form of Hf'^+ ion, 
similar to the Zr^^ centres observed in zirconia'*'*. The 
electron self-trapping if effective, may present a serious 
problem in the high-k MOSFETs. However, our density 
functional calculations do not predict this defect as sta- 
ble, possibly due to the known limitations of local approx- 
imations regarding localisation problem^. Therefore, we 
leave the problem of small polaron assisted trapping for 



further investigation. 

C. Fixed charge effects 

Both N and N2 intcrstitials and subsitions are deep 
centres, and according to Figure |S1 their relaxed affinities 
(electrical levels) are all below the VBM of Si, which 
would make them type C defects. This means that 
at equilibrium they ought to have maximum negative 
charge. However, such an equilibrium state assumes that 
the electrons are supplied from the silicon valence band, 
which is difficult (if not impossible) for two reasons: 1. 
The effectiveness of the resonant tunneling strongly de- 
pends on the overlap between the states involved in the 
process. Therefore, defects localized far away from the 
interface are ineffective electron traps, even if their en- 
ergy is resonant with Si VB. 2. Charging occurs predom- 
inantly via resonant electron tunneling from the Si VB, 
i.e. it envolves transitions on to the unoccupied defect 
states. These states (shown as unrelaxed levels in Figure 
|HJ| are mostly of type B. Therefore, most of the nitrogen 
species are unavailable for a resonant tunneling from the 
silicon valence band. Hence, they may retain multiple 
non-equilibrium charge states for a very long time, and 
thus contribute to either positive or negative charge. The 
possible exceptions are NO intcrstitials and N substitu- 
tions, whose electrical levels are below the Si VBM, so 
they will quickly achieve maximum negative charge. 

The only definite source of positive fixed charge in our 
calculations are protons and anion vacancies. With re- 
gard to these defects, atomic nitrogen and molecular ni- 
trogen will reduce vacancy concentration, but may them- 
selves produce negative or positive charge states as diss- 
cussed above. Ammonia anneal will result in higher con- 
centration of protons, but it is unclear how much of these 
will be evacuated during the anneal, so the increase in 
positive charge is expected with the NH3 PDA. 

Finally, we would like to emphasise that high-k mate- 
rials generally favour charged defect states, due to their 
large polarization energy. This is quantified for m-Hf02 
in Table in 

We conclude that nitrogen's capacity for reducing fixed 
charge in hafnium oxide is limited. Therefore, the prob- 
lem of fixed charge built up and Vt instability is not 
resolved by PDA. Hinkle and Lukowsky^^ reported that 
bulk only nitridation does not decrease and may even 
increase fixed charge in the oxide film, which is in quali- 
tative agreement with our calculations. 

D. Growth of the interfacial SiOa; layer 

We assume that growth of the SiO^, layer at the inter- 
face is controlled by a migration of oxygen towards the 
interface. Therefore, we considered how various PDA 
gases interact with oxygen intcrstitials. Our calculations 
predict N2 does not bond to oxygen interstitial in any 
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FIG. 8: Electron affinities of various defects in hafnia, given in botli vertical and deeper, relaxed forms. All values are in eV. 



charge state (Table ^Oj. Therefore, interstitial oxygen is 
free to diffuse towards the interface. In contrast, atomic 
N strongly bonds to the interstitial oxygen, forming less 
mobile NO^ molecules. As a result oxygen diffusion to- 
wards the interface is blocked, thus inhibiting growth of 
interfacial SiO^: layer. This conclusion is consistent with 
data on nitrogen remote plasma deposition processing!^, 
where reduction of interfacial silica growth has been ob- 
served. 

Ammonia molecules in m-Hf02 dissociate sponta- 
neously resulting in the formation of NHJ^~, NH*^, H+ 
and possibly H2 and H^. Of these products, only pro- 
tons will effectively bond to the interstitial oxygen ions 
and exclusively to those in the charge state -2. At the 
same time, the most mobile charge state of oxygen, 0~— , 
does not interact strongly with protons. Therefore, am- 
monia PDA will not inhibit effectively the growth of the 
interfacial SiOa; layer. The same applies to the PDA with 
nitric oxide. 



ited effect on the channel mobility. So, to understand 
the processes in the channel, an extension of this work is 
required towards systems which include interface effects 
- work in this direction is currently under way. Also, 
we note that bulk DFT calculations only give an indi- 
cation of the defect charge states that are thermody- 
namically possible. Some of these states may or may 
not be relevant depending on the values of the electron 
chemical potential. Moreover, charge transport through 
the dielectric is generally slow, so actual defect states 
may be controlled kinetically rather than thermodynam- 
ically. Additionally, calculations of relative energies of 
defects in different charge states require empirical cor- 
rections due to various limitations of the method applied 
(e.g ill-defined zero energy and spurious electrostatic in- 
teractions in periodic boundary conditions, underestima- 
tion of the energy band gap etc.). Such corrections must 
be treated with some caution, and new theoretical devel- 
opments are required for a more rigorous approach. 



E. Perspective and future work 

Our calculations suggest that hafnium oxynitride 
might have better properties as a gate oxide. Due to 
the multiplicity of the NO oxidation states, this mate- 
rial may form stable amorphous structures with accept- 
able dielectric properties. Some preliminary results on 
hafnium oxynitride gates in MOSFET have been recently 
reported^, but a more systematic study of film com- 
position and electrical properties is required. 

In this work we have considered only bulk structures 
close to their thermodynamic equilibrium. There exists 
strong experimental evidenc o^^'^^ that in certain cases a 
significant part of the electron traps are located in the 
bulk of the insulator. However, fixed charge distant from 
the channel is well screened and therefore, has only lim- 
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APPENDIX A: CALCULATION OF FORMATION 
AND IONIZATION ENERGIES, AND ELECTRON 
AND HOLE AFFINITIES 

The formation energy of a defect in the charge state q 
is given by: 

G{D\Ef) = El- El - E'j^' + qEp, (Al) 

where E'^ denotes the total energy calculated for the sys- 
tem X at the equilibrium geometry and the excess super- 
cell charge q in units of the electron charge e, E^ is the 
energy of the perfect Hf02 crystal in the charge state g, 
E^j^"^ is the energy of the isolated molecule, and all en- 
ergies are calculated in the same supercell. Depending 
on the value of Ep, the thermodynamically stable charge 
state of the specific defect may change. 

In order to study stable charged defect states and the 
possible role of defects in photo- and thermo-stimulated 
processes, as well as in electronic devices one needs to 
know the electron affinities and ionization energies of the 
defect states with respect to the bottom of the conduc- 
tion band of hafnia and to other electron or hole sources, 
such as silicon. To achieve that, we compare total energies 
of the initial and final systems with the same number of 
electrons. The main inaccuracy of this approach is re- 
lated to the underestimation of the band-gap in the DFT 
calculations. This means that the defect states appear to 
be closer to the valence band and conduction band edges 
in hafnia. The relative error in the energy level position 
with respect to the gap edges depends on a defect and is 
impossible to establish without proper calibration using 
experimental data. 

Defining the absolute value of the defect ionization en- 
ergy Ip{D'^) as the vertical excitation energy of an elec- 
tron from the defect with the charge q to the bottom of 
the conduction band, we have: 

/p (D'^ ) = E^ +El+'-E",-El + Kl, ( A2) 

where the value E^^^ is calculated for the geometry of 
the relaxed defect with charge q and k1 is a correction 
for the position of the bottom of the conduction band. 
Similarly we can define the electron affinity of the defect 
Xe{D'^) (i-e. the energy gain when the electron from the 
bottom of the conduction band is trapped at the defect) 
as follows: 

XeiD") ^ E^ + El ~ E''^ - E^,^^ + k2. (A3) 

Here the correction k2 can be generally different from k1. 
One can consider both "vertical" and "relaxed" electron 



affinities. In the latter case the lattice relaxation after 
the electron trapping is included in El^^ . We can also 
define the hole affinity of the defect Xh{D'^), i-e. the 
energy gain when the a free hole is trapped from the top 
of the valence band to the defect as follows: 

Xh {D^ )^E++El^ El -E%+'+k3. ( A4) 

Again, dependent on whether the lattice relaxation in 
the final state is included or not, one will obtain different 
affinities. The vertical hole affinity provides a useful es- 
timate of the position of the defect state with respect to 
the top of the valence band. One can also verify, that the 
relaxed affinity is related to the defect's q-|-l/q electric 
level as defined elsewhere (e.g. Ref42i). 

To define the corrections nl, k2, k3 we use the fol- 
lowing considerations, i) We assume that the main inac- 
curacy in defining the relative positions of defect states 
with respect to the band-gap edges is due to unoccupied 
Kohn-Sham states, and that the underestimated band 
gap is mainly due to the too low position of the conduc- 
tion band minimum. Therefore, we use an approximation 
that k1 = k2 — K and k3 = 0. These conditions are diffi- 
cult to fully justify without comparison with experiment, 
ii) Using these conditions and definitions HA3|I and (|A4(I 
it is easy to obtain: 

Xh{D'^) + X.{D^^')^E,{exp), (A5) 

where both affinities correspond to relaxed final defect 
states. This condition holds in all calculations, which 
insures the consistency of our approach, iii) We use the 
experimental value of Eg{exp) = 5.68 eV— to define the 
difference 



K = Eg{exp) — Eg{theor), (A6) 

and correct the defect ionization potentials and electron 
affinities. This gives k = 5.68 — 3.92 = 1.76 eV, which is 
used in all calculations. 

Although this method is approximate, fixing the value 
of K allows us to present the results of our calculations 
in one scale. Another advantage is that, in order to find 
defect affinities with respect to electrons at the bottom 
of silicon conduction band or holes at the top of silicon 
valence band, within the same method one can use the 
experimental value of band offset with Si. This scale can 
be changed if a more "accurate" or relevant value for k will 
be found. This will require only a shift of our predicted 
values by a constant. 
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